Stabilization of in-phase fluxon state by geometrical confinement in small 

Bi2Sr2CaCu208+^ mesa structures. 
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The in-phase (rectangular) fluxon lattice is required for achieving coherent THz emission from 
stacked Josephson junctions. Unfortunately, it is usually unstable due to mutual repulsion of fluxons 
in neighbor junctions, which favors the out-of-phase (triangular) lattice. Ifere we experimentally 
study magnetic field modulation of the critical current in small Bi-2212 mesa structures with different 
sizes. A clear Fraunhofer-like modulation is observed when the field is aligned strictly parallel to 
superconducting CuO planes. For long mesas the periodicity of modulation is equal to half the flux 
quantum per intrinsic Josephson junction, corresponding to the triangular fluxon lattice. However, 
the periodicity is changed to one flux quantum, characteristic to the rectangular fluxon lattice, both 
by decreasing the length of the mesas and by increasing magnetic field. Thus, we demonstrate that 
the stationary in-phase fiuxon state can be effectively stabilized by geometrical confinement in small 
Bi-2212 mesa structures. 

PACS numbers: 74.72.Hs, 74.78.Fk, 74.50.+r, 85.25.Cp 



Josephson flux-flow oscillators can provide a remark- 
able linewidth ^ 1 Hz in the sub-THz frequency range 
but with a small power < fiW . The emission power 
can be increased by phase-locking of several coupled os- 
cillators d, The strongest coupling is achieved be- 
tween atomic scale intrinsic Josephson junctions (IJJs), 
naturally formed in single crystals of high temperature 
superconductor Bi2Sr2CaCu2 08-i-x (Bi-2212) [4|. Fur- 
thermore, a large energy gap in Bi-2212 would al- 
low operation in an important THz frequency range. 
Therefore, stacked IJJs are intensively studied as pos- 
sible candidates for tunable, high-power THz oscillators 

aii0,ii[ioi. 

To achieve power amplification, junctions must be 
locked in the in-phase mode, so that Josephson vortices 
(fluxons) form a rectangular lattice in the stack. Unfor- 
tunately, the rectangular lattice is usually unstable due 
to mutual repulsion of fluxons, which favors formation 
of the triangular (out-of-phase) fluxon lattice. Motion of 
the triangular lattice leads to out-of-phase oscillations in 
neighbor junctions, which results in destructive interfer- 
ence and negligible emission. Thus, the major challenge 
for achieving coherent THz emission from stacked IJJ is 
to promote the in-phase fluxon state. 

Fluxon distribution in stacked junctions is governed 
by three forces: (i) the in-plane fluxon-fluxon repulsion; 
(ii) the in-plane fluxon-cdge interaction and (iii) the in- 
terlayer fluxon-fluxon repulsion. The latter is specific 
for stacks and promotes the triangular lattice. The two 
in-plane forces promote the rectangular fluxon lattice. 
Therefore, the rectangular lattice can be stabilized via 
decreasing the interlayer force by decreasing the inter- 
layer coupling PH . or via increasing the in-plane forces 
by reduction the junction length L [1,0, [U or by increas- 



ing the fluxon density [1, 0| ■ The latter can be described 
by the approximate inequality: 



$/$o > L/\j 



(1) 
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where \,j the Josephson penetration depth, $ = HLs is 
the flux through each IJJ, s ~ 1.5nm is the interlayer 
spacing of Bi-2212 and $o = hc/2e is the flux quantum. 

The symmetry of fluxon lattice can be revealed from 
modulation of the critical current as a function of in- 
plane magnetic field Ic{H) [7]: the rectangular lattice 
leads to conventional Fraunhofer Ic{^) modulation with 
$0 periodicity, while the triangular lattice leads to 
periodicity because the fluxon lattice parameter in the 
c— axis direction is doubled. Previous measurements of 
Ic(H) for IJJs didn't reveal clear periodic modulation 
[1 [ill, EjJEjJE- The same is true even for low-T^ 
stacks [j, [l^ [201 ■ This was explained by existence of 
multiple metastable fluxon modes in long, strongly cou- 
pled stacked junctions [ii[,[il|. So far, clear $o/2 and $0 
modulations were observed onhMn the dynamic flux-flow 
state both for low-T^ stacks 15] and IJJs [H, [H, [H. 
Yet, their interpretation is ambiguous because even the 
rectangular lattice may lead to the $0/2 modulation in 
the flux-flow state [l^. Therefore, it is necessary to ana- 
lyze the static critical current, not the least because this 
is the most important fingerprint of the dc- Josephson ef- 
fect. 

Here we study size and magnetic field dependence of 
the critical current in small Bi-2212 mesa structures. We 
observe a clear Fraunhofer-like modulation of Ic as a 
function of the in-plane magnetic field. The periodic- 
ity of modulation changes from half $0 - characteristic 
for the triangular lattice, to $0 - characteristic for the 
rectangular lattice, both upon increasing the field and 
decreasing the junction size. This provides a clear exper- 
imental evidence that the in-phase fiuxon state can be 
stabilized by geometrical confinement in small stacked 
Josephson junctions, even in the stationary case. 



Mesa 


L [/im] 


W [^m] 


Jco [mA] 


Aj [/im] 


J^o [T] 


i/o""'' [T] 


4a 


5.1 


1.4 


75.0 


0.68 


0.27 


0.26 


4b 


2.7 


1.4 


39.5 


0.69 


0.51 


0.55 


6b 


2.0 


1.7 


39.8 


0.65 


0.69 


0.662 



TABLE I: Properties of the measured mesas on the same Bi- 
2212 single crystal. W is the width of the junctions parallel 
to the magnetic field, Aj = {^ocsLW) /{IQn'^IcO^'ii,), with 
s — 1.5 nm and Xat ~ 200 nm and Ho ~ ^o/Ls, H^^^ is the 
measured periodicity of /c-oscillations. 



Mesa structures, containing few IJJs, were fabricated 
on top of Bi-2212 single crystals with Tc = 82K. Some 
mesas were FIB-trimmed to reduce the size. Details of 
mesa fabrication can be found in Ref. [2l|. Table-U 
summarizes properties of studied mesas. Samples were 
mounted on a rotatable sample holder with alignment ac- 
curacy ^ 0.02°. All measurements reported below were 
performed at T — 1.6 K. 

In Bi-2212 mesas the top surface junction is always de- 
teriorated and has much smaller Ic than the rest "bulk" 
IJJs [l^]. With increasing current, /, it switches into the 
quasiparticle branch Vi(/), while "bulk" IJJs remain in 
the stationary superconducting state up to much larger 
"bulk" critical current. Since we are interested in the 
collective behavior of the stack, associated with the sym- 
metry of the fluxon lattice in the stack, we have to study 
the "bulk" critical current. To do so, we first carefully 
measured the quasiparticle branch of the surface junc- 
tion Vi(/) at iJ = 0. We then automatically subtracted 
it from the I — V characteristics of the mesas. 

Fig. [T] shows current voltage characteristics of mesa 
4b with the subtracted first quasiparticle branch Vi, for 
different magnetic fields. It is seen that subtraction works 
very well, with accuracy better than lOfiV at all fields. 
The bulk critical current Ic was determined from the first 
deviation of V — Vi from zero, using a 100/iV criterion. 
A strong modulation of the critical current as a function 
of magnetic field is evident from Fig. [TJ With increa sing 
field, the flux-flow branch with distinct Fiske steps |14| 
develops in the I — V curve, as seen from Fig. [T]b). 

Field alignment was crucial for observation of Ic{H) 
modulation. Even slight misalignment with respect to 
CuO planes led to avalanche entrance of Abrikosov vor- 
tices with increasing field. After that the Ic{H) pat- 
terns become heavily distorted, irreversible and com- 
pletely lack periodic modulation, as shown in inset of 
Fig. O Thus, exact field alignment was the main exper- 
imental challenge in this work. 

At low fields H < 6T the alignment was complicated 
by field lock-in [Ts] along the CuO-plane, which led to a 
hysteresis of as much as three degrees, making the precise 
alignment impossible. To avoid the hysteresis, we looked 
at the high field H — 15T magnetoresistance at bias close 
to the sum-gap voltage [13]. The latter is very sensitive 
to the c-axis field component, which leads to a substan- 
tial negative magnetoresistance [1^. To the contrary. 




FIG. 1: I — V characteristics of mesa 4b with the subtracted 
quasiparticle branch of the surface junction, at different mag- 
netic fields: a) within the central Meissner lobe of /c(*&) and 
b) for 1 < $/$o < 2. Flux-flow and Fiske steps are seen. 



the magnetoresistance was negligible for precise in-plane 
field alignment. Thus we were able to accurately align 
samples by minimizing the high field and high bias mag- 
netoresistance. Even after perfect alignment, presence of 
trapped Abrikosov vortices affected the Ic, as shown by 
red line in Fig. 2. Yet, immobile Abrikosov vortices did 
not distort the qualitative shape of Ic{H) modulation. 
To get rid of trapped Abrikosov vortices, samples were 
always heated up above Tc after alignment. 

Figs. [2] and [3] present the main experimental result of 
our work: magnetic field modulation of the "bulk" criti- 
cal current for Bi-2212 mesas of different sizes. Vertical 
grid lines correspond to integer flux quanta per I J J. All 
shown patterns are perfectly reversible. Magnetic field 
was swept back-and-forth several times (from -2 to +2 
T in Figs. [Hand Ob), and from -4T to 4T in Fig. [3] a) 
without significant changes in IdH). This confirms that 
alignment was fine and Abrikosov vortices didn't enter 
the mesa up to the largest field. 

Figure [5] a) shows the Ic{^) pattern (in a semi-log 
scale) for the long mesa (4a), with L = 5.1 /im and 
L/Xj ~ 7.5. When increasing the field from zero, Ic first 
continuously decreasing and a,t H > 0.5 T, the oscillatory 
behavior appears. The period of modulation is initially 
$0/2 rather than $0j with maxima both at integer and 
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FIG. 2: Experimental Ic(^) for the long mesa: without 
Abrikosov vortices (blue line, four field sweeps) and with 
few immobile Abrikosov vortices (red line, two field sweeps). 
Three field regions can be distinguished: i) "I> < 2<l?o: 
metastable region with strong fluctuations and without mod- 
ulation; ii) 2$o < $ < 5$o: 'I'o/2 modulation with even in- 
teger and half- integer maxima; iii) $ > 5<l?o: transition to "l>o 
modulation with predominance of half-integer maxima. Inset 
shows /c("I>) for two sweeps down from lOT. 
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half-integer <I>o- However, with increasing field the am- 
plitude of integer maxima rapidly decreases and at high 
fields maxima at half-integer <i>/$o become dominant. 

Figure [3] a) shows Ic{^) (in a semi- log scale) for the 
moderately long mesa 4b with L = 2.7 ^uy and L/Xj ~ 
3.9. Here maxima at integer $/$o are strongly sup- 
pressed even at low fields. At the condition of Eq.iJT]), 

> 3$0j the integer maxima completely disappear and 
the Ic{^) switches to conventional Fraunhofer modula- 
tion with one flux-quantum periodicity and maxima at 
half-integer ^/^o- 

Finally, /c(*&) for the shortest mesa (6b) is shown in 
figure [3] b). The length of the mesa is L = 2.0 /im and 
L/Xj = 3.1. Here integer maxima are suppressed almost 
completely already for $ = 2<i>o. 

For better understanding of experimental data, in Fig. 
|3]we show numerically simulated Ic{^) for N — 5 stacked 
IJJs of different lengths: a.) L ~ 8Xj (in the semi-log 
scale), b) L = AXj and c) L — 2Xj (for details see the 
Supplementary [2J])- It is seen that the overall behavior 
of Ic{^,L) is very similar to experimental data in Figs. 
[2] and [2 

For the long stack in Fig. 4 a) three field regions can 
be distinguished, similar to that in Fig. 2: 

i) At low field <f> < 2$o there are strong fluctuations of 
Ic due to chaotic switching between multiple metastable 
fluxon modes [HiliBl, and no clear modulation of /c($). 

ii) At intermediate fields 2$o < $ < 6^o modulation 
with periodicity of <I>o/2 appears. The fluxon lattice is 
not yet formed in this fleld range and modulation is due 



FIG. 3: Experimental Ic{^) for a) moderately long (in a semi- 
log scale) and b) the shortest mesas. It is seen that integer 
maxima are rapidly decreasing with increasing field and de- 
creasing L. At ^/$o > L/Xj only half-integer maxima are 
seen, indicating transition to the rectangular fluxon lattice. 
Inset in a) shows metastable sub-branches at "I>/$o ~ 1- 



to switching between certain quasi-periodic modes '24|. 

iii) Finally at high fleld, approximately determined by 
condition of Eq. ([T]) , the regular fiuxon lattice is formed 
and maxima at integer <&/$0j corresponding to the trian- 
gular lattice become weaker than maxima at half-integer 
^/^O: corresponding to the rectangular lattice, so that 
periodicity of modulation becomes $0. 

In shorter stacks, the metastable region becomes 
smaller. In Figs. 3 a) and 4 b) it is seen only at 
$ ^ $0- In inset of Fig. 3 a) we show this metastable 
region, measured by repeated field sweeps. Certain sub- 
branches can be distinguished, indicating a reduced vari- 
ety of metastable fluxon modes in comparison to the long 
junctions case. In Figs. 4 b) we marked the most prob- 
able modes, which cause distinct subbranches in Id^)- 
The metastable states gradually die out with decreasing 
L, as seen from the data for the shortest junctions in 
Figs. 3 b) and 4 c), and eventually disappear for L < Xj. 

Even region-ii) disappears in shorter stacks and the 
amplitude of integer maxima is always much smaller than 
of half-integer maxima. The integer maxima are further 
suppressed with increasing field, see Figs. 4 b) and 3 
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FIG. 4: Simulated Jc(<E>) for stack of iV = 5 IJJs for a) 
L = 8Aj (in a semi-fog scale), b) L = 4Aj and c) L = 2A,j. 
The transition from $o/2 to "l?o periodicity occurs with in- 
creasing <1> and decreasing L and reflects the transition from 
triangular to rectangular fiuxon lattice, as shown by snap- 
shots in panel a). Shaded areas indicate regions of stability 
of the rectangular (in-phase) lattice. The in-phase mode is 
dominating in the shortest stack so that Ic{^) returns to the 
conventional Fraunhofer pattern, shown by the blue dashed- 
dotted line. 



a) and with decreasing junction length, see Figs. 4 c) 
and 3 b). For shorter junctions the $o periodicity is es- 
tabhshed at smaller $/$o, in agreement with Eq.(|T]). In 
short stacks, L < Xj, only maxima at half-integer $/<i>o 
are left, and Ic{^) returns to the conventional Fraun- 
hofer pattern, shown by the blue dashed-dotted line in 
Fig. 4 c). The apparent similarity between measured 
and simulated data in Figs. 2-3 confirms our estimation 
of A,7 - 0.7 fim in Table- |I| 

The observed transition from <f>o/2 to (f>o periodicity of 
/c(<&) modulation is the consequence of transition from 
triangular to rectangular fiuxon lattice in the stack 0] ■ 
This is illustrated by snapshots in Fig. 4 a), which show 
fiuxon distributions at both sides of the subdominant 
maximum at $/$o — 7 and the dominant maximum at 
$/$o — 7.5. It is seen that the rectangular (in-phase) lat- 
tice occurs at the dominant half-integer maximum. For 
the stack with L = 2Aj the in-phase state is stable prac- 
tically in the whole field range, marked by shaded area 
in Fig. 4. All this could be verified using the real-time 
simulation code provided in the Supplementary . 

In summary, we observed clear Fraunhofer-like mod- 
ulation of Ic{H) in Bi-2212 mesa of different sizes. We 
demonstrated that the periodicity of modulation changes 
from half-flux quantum to flux quantum both with in- 
creasing magnetic field and decreasing mesa size. This 
indicates transition from triangular (out-of-phase) to 
rectangular (in-phase) fiuxon lattice in stacked intrinsic 
Josephson junctions. Therefore, we confirm experimen- 
tally that the stationary in-phase fiuxon state can be 
stabilized by geometrical confinement in small Bi-2212 
mesa structures. This is important for realization of high 
power fiux-fiow oscillator in the THz frequency range. 

We are grateful to A.Tkalecz and A.Rydh for assis- 
tance in experiment. The work was supported by K. & 
A. Wallenberg foundation, the Swedish Research Council 
and the SU-Core Facility in Nanotechnology. 



[1] V.P.Koshelets and S.V. Shitov, Supercond. Sci.Technol. 

13, R53 (2000). 
[2] P.Barbara, et al., Phys.Rev.Lett. 82, 1963 (1999) 
[3] S.Sakai, et al., Phys.Rev.B 58, 5777 (1998) 
[4] R.Kleiner et al., Phys.Rev.B 50, 3942 (1994) 
[5] N.Gr0nbech- Jensen, et al., Phys.Rev.B. 50, 6352 (1994) 
[6] M.Machida, et al., Physica C 330, 85 (2000); M. Machida, 

Phys. Rev. Lett. 96, 097002 (2006). 
[7] A. E. Koshelev, Phys. Rev. B 66, 224514 (2002); ibid 75, 

214513 (2007). 
[8] I.E.Batov, et al., Appl.Phys.Lett. 88, 262504 (2006) 
[9] L. Ozyuzer, et al., Science 318, 1291 (2007). 
[10] M.H.Bae, et al., Phys.Rev.Lett. 98, 027002 (2007) 
[11] V. M. Krasnov et al., Physica C 304, 172 (1998). 
[12] V.M.Krasnov, et al., Phys.Rev.B 54, 15448 (1996) 
[13] S. Ooi, T. Mochiku, and K. Hirata, Phys. Rev. Lett. 89, 

247002 (2002). 
[14] S.M.Kim, et al, Phys.Rev.B 72, 140504(R) (2005) 



[15] K.Kadowaki, et al., Physfoa C 437-438, 111 (2006) 
[16] A.V.Ustinov and N.F.Pedersen, Phys.Rev.B 72, 052502 
(2005) 

[17] Y.I. Latyshev, J.E. Nevelskaya, and P. Monceau, Phys. 

Rev. Lett. 77, 932 (1996). 
[18] V.M.Krasnov, et al., Phys.Rev.B 61, 766 (2000) 
[19] "l>o modulation of Jc was reported in unconventional "hol- 
low" junctions by Y.I. Latyshev et al, Physica C 362, 156 
(2001); the $o/2 modulation similar to that in Fig. 2 was 
reported by I. Kakeya, et al., Physica C 437-438, 118 
(2006) 

[20] H.Kohlstedt, et al., IEEE Trans. Appl. Supercond. 3, 
2197 (1993); H.Amin, M.G.Blamire and J.E.Evetts, ibid, 
3, 2204 (1993). 

[21] V.M.Krasnov, T.Bauch and P.Delsing, Phys.Rev.B 72, 

012512 (2005) 
[22] V.M.Krasnov, Cond-Mat 0902.3588 

[23] N. Morozov et al., Phys. Rev. Lett. 84, 1784 (2000); V.M. 



5 



Krasnov, ot al., ibid 86, 2657 (2001). the simulation program available in the Supplementary: 

[24] The variety of fluxon states can be explicitly seen form see EPAPS Document No... 



